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Abstract 

^ ■ Measured hadron yields from relativistic nuclear collisions can be equally 



well understood in two physically distinct models, namely a static thermal 



^ ■ hadronic source vs. a time-dependent, nonequilibrium hadronization off a 



quark-gluon plasma droplet. Due to the time-dependent particle evapora- 
tion off the hadronic surface in the latter approach the hadron ratios change 
(by factors of ^ 5) in time. Final particle yields reflect time averages over 
the actual thermodynamic properties of the system at a certain stage of the 
evolution. Calculated hadron, strangelet and (anti-) cluster yields as well as 
freeze-out times are presented for different systems. Due to strangeness dis- 
tillation the system moves rapidly out of the T, plane into the /Hs-sector. 
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Strangeness to baryon ratios /s = 1 — 2 prevail during a considerable frac- 
tion (50%) of the time evolution (i.e. A-droplets or even S~-droplets form 
the system at the late stage: The possibility of observing this time evolu- 
tion via HBT correlations is discussed). The observed hadron ratios require 
Tc ^ 160 MeV and B^/^ ^ 200 MeV. If the present model is fit to the extrap- 
olated hadron yields, metastable hypermatter can only be produced with a 

probability p < 10~^ for A>4. 
25.75.Dw, 12.38.Mh, 24.85.+p 
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I. HADRON RATIOS FROM AN EQUILIBRATED SOURCE 



Measurements of strange and anti-strange particles in relativistic nuclear collisions have 
received much attention recently [0,0: (Anti-) strangeness enhancement relative to pp-data 
has been predicted as quark gluon plasma signal, because equilibrium yields of strange 
hadrons cannot be achieved during the short collision time due to small ss cross sections 
0]. However, several thermal models have been apphed [^-^ to (strange) particle yields and 
used to extract the characteristic thermodynamic properties of the system (a few macroscopic 
parameters) from chemical equilibrium. 

Fig. ID shows such a calculation: different measured and calculated hadron ratios are 
shown as resulting from an equilibrium hadron gas at fixed temperature T, quark chemical 
potential fig and strange quark chemical potential fis- The constituents of the hadron gas 
are all well-established hadrons up to masses of 2 GeV (coinciding with the list given in [0) 
taken from [§. Note that an isospin symmetric system is assumed (as in the hadronization 
model): # (protons) =#(nucleons)/2, #(7r+)=#(all pions)/3 etc. The system is treated as 
a mixture of relativistic, non-interacting Bose-Einstein and Fermi-Dirac gases with Hage- 
dorn's eigenvolume correction. 

The parameters T and fig of Fig. |1| are chosen in line with the values given by P. Braun- 
Munzinger and J. Stachel in |Q. Agreement with many observed ratios from different ex- 
periments at the SPS is achieved (as in ^j). The value of /i^ has been determined by 
the requirement that the total net-strangeness vanishes: = 0. The quantity fs is the 
fraction of net-strangeness over net-baryon number {fs = {Ng — Ns)/A). Note that the 
presently calculated strange quark chemical potential fis = 24 MeV is different from the 
value /is = 18.6 MeV given in 0. For the ratios we gain roughly the same results, however, 
there is a discrepancy of about a factor of 2 in the A/{p — p) and A/p values. 

It should be mentioned that some of the measured ratios may deviate from the true, yet 
to be observed 47r-ratios: the p/ir^ ratio, which is a rather direct measure of the specific 
entropy S/A in the fireball, has a value of 0.18(3). However, it is calculated from proton 
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yields in the rapidity range 2.6 to 2.8 and pion yields in the range 3.2 to 4.8 0. The 
+ Q + S ) ratio is measured for pt > 1.6 GeV which excludes a major part of 

the momentum distributions. If one takes the value for a common mj-cut (1.7 ± 0.9 for 
rrit > 2.3 GeV ||T^) the ratio is compatible with the extracted temperatures and chemical 
potentials, since within the restricted phase space it is just the ratio of the fugacities times 
the degeneracies [^], which deviates considerably from the value for the full phase space. 
However, to allow for a consistent comparison to other models all experimental data are 
taken from the compilation in ^ and we also give Air ratios for both presented models. 
Fig. ID shows the massive effect (factors of ^ 3!) of feeding on the predicted ratios (cf. 



the discussion in ||TTI]): the crosses denote the ratios as predicted, if the contributions due to 
the decay of higher lying resonances (e.g. A ^ n + p, p + etc.) are ignored. The 

circles show the ratios as obtained after one includes these feeding effects. In some cases 
both hadron species change considerably due to feeding, but sometimes such that the ratio 



does not change much, e.g. {p/pY'^ ~ ^eg _|_ , where p*, p* denote the yields [|1T| due to 
the decay of resonances (e.g. A, A). 

Collective motion and fits to particle spectra, respectively, are not considered here. Lon- 
gitudinal and transverse flow velocities have already been extracted from particle spectra in 
combination with a static thermal model |]12|. However, it is questionable, whether thermal 
and chemical freeze-out happen simultaneously. Here, we focus on a collective, chemical 
freeze-out. 



II. THE HADRONIZATION MODEL 

We now investigate, whether the reasonable agreement of data and model, as reached 
with the static equilibrium-plus-feeding model, can also be achieved with a dynamic model, 
which includes the formation and expansion of quark-gluon matter, a first order phase 
transition into coexistence of quark and hadron matter, and time-dependent evaporation of 
hadrons from the system, as it evolves with time through the phase transition. This means 
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adopting a model for the hadronization and space-time evolution of quark matter droplets 
given in ||13|. A similar approach was given in [Q, however, in terms of rate equations for 
the flavor kinetics. They allow the abundances of the different constituents of the system 
to be out of equilibrium, which is in contrast to our model. The calculated hadron yields in 
both models are similar, if the same initial input is used. 

The model assumes a first order phase transition of the QGP to hadron gas with 
Gibbs conditions (P<3GP = pHG ^ j.qgp ^ j.hg ^ ^qgp ^ ^hg ^ ^qgp ^ ^hg^ coexistence. 

The expansion and evaporation of the system takes into account equihbrium, as well as 
nonequilibrium features of the process: 

1. The plasma sphere is permanently surrounded by a layer of hadron gas, with which it 
stays in thermal and chemical equilibrium during the phase transition (Gibbs condi- 
tions). The strangeness degree of freedom stays in chemical equilibrium between the 
two phases (however, < s — s >^ for the individual phases). Thus, the hadronic 
particle production is driven by the chemical potentials. 

2. The particle evaporation is incorporated by a time- dependent freeze-out of hadrons 
from the hadron phase surrounding the QGP droplet. During the expansion, the 
volume increase of the system thus competes with the decrease due to the evaporative, 
time-dependent freeze-out. 

The global properties, like (decreasing or increasing) S/A and fs of the remaining two- 
phase system, then change in time according to the following differential equations for the 



baryon number, the entropy, and the net strangeness number of the total system |jT3| 

d 



^\tot -p aHG 

d 

^(iv,-iv^)*°* = -r3(iv,-iv^)^^ 



dt 



For simplicity, we set Fi = r2 = F3, thus, F = -^^u ( ^^t° ) effective ('universal') 

rate of particles (of converted hadron gas volume) evaporated from the hadron phase. A 



more general treatment with differing rates for the three quantities is presently being studied 

The equation of state consists of the bag model for the quark gluon plasma and the 
hadron gas of the previous section. Thus, one solves simultaneously the 'non-equilibrium' 
dynamics (|I]) and the Gibbs phase equilibrium conditions for the intrinsic variables, i.e. the 
chemical potentials and the temperature, as functions of time. 



III. PARTICLE RATES FROM THE HADRONIZING PLASMA 

The particle yields as functions of time have been calculated for different parameter 
sets: In Fig. ^ the final (time integrated) ratios are plotted for a bag constant of B^^'^ = 
235 MeV, initial strangeness fraction = and an initial specific entropy per baryon of 
S/A = 45. We choose an initial net-baryon number A^** = 100. Obviously, the particle 
ratios will not depend on this choice. The theoretical ratios show an equally good overall 
agreement with the data points as the static fit. This scenario results in a rather rapid 
hadronization. The quasi- isentropic expansion of the system is due to those hadrons which 
are subsequently evaporated. Decays and feeding occur after the hadrons leave the system. 
The specific entropy extracted with this fit is higher than the values given in for this 



beam energy {S/A ^ 20 — 25), where entropy production was calculated in a one-dimensional 
hydrodynamic model of a heavy ion collision. The calculation of particle rates from a 
hadronizing QGP with B^^'^ = 235 MeV has already been presented in [1^, however, with 
S/yl*™* = 25. Predictions have been also made for different conditions within the rate- 



equation approach of [|Tj] which yields much higher K/ir ratios. 

The system can cool within this model (meaning the employed equations of state for 
both phases) if {S/A)'^'^^/{S/A)'^^'^ < 1 [|1^,|1^], which can only be achieved for a low bag 
constant. Low bag constants allow for the formation of metastable strange remnants of the 
plasma, the strangelets. However, the particle yields as resulting from such small -B-values 
in the present model are not in good overall agreement with the observed ratios. This is 
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shown in Fig. |^ where a bag constant of 5^/^ = 160 MeV and an initial specific entropy per 
baryon of S/A = 40 and S/A = 150, respectively, is chosen. The results for the lower specific 
entropy are tolerable except for ratios where antibaryon yields are set in relation to baryon 
or meson yields. Antibaryon to baryon ratios fail by more than one order of magnitude, 
which is due to the lower temperatures of the coexistence phase. To counteract this effect, 
much higher initial specific entropies would be required. However, in particular the ratios 
p/vr"*" and d/p will then move away considerably from the experimentally observed values 
(see Fig. m. 



Note that within microscopic hadronic model calculations |T8[ a pronounced dependence 
of antibaryon (and even more so anti-cluster) yields on the reaction volume has been pre- 
dicted. This chemical non-equilibrium refiects the strong sensitivity of antibaryon produc- 
tion and absorption on the phase space evolution of the baryons. It is in contradiction to 
the 'volume freeze-out' of thermal models. On the other hand, even if equilibration of an- 
tibaryons is assumed, the abundances are predicted to be strongly enhanced in the case of 



an interacting hadron gas [|T9 



Fig. ^ shows the time evolution of the system (initial conditions as in Fig. |^) in the 
plane of quark and strange quark chemical potential. The quark chemical potential is at the 
beginning of the evolution in the order of the temperature of the system. It is /ig ~ 110 MeV. 
This is twice the value assumed in the static approach, fig ^ 60 MeV. However, in the course 
of the hadronization, the quark chemical potential drops to ~ 15 MeV. On the other hand, 
the strange chemical potential increases from /i^ = to values of fig ~ 50 MeV at the end of 
the evolution, as compared to /is ~ 25 MeV in the static fit. This is due to the strangeness 
distillation effect (see below). 

Fig. 1^ shows the time evolution of the system in the /ig-T plane in conjunction with the 
phase diagram of the quark gluon plasma phase and the hadron gas. The initial specific 
entropy oi S/A = 45 in combination with the bag constant 5^/^ = 235 MeV leads to a fast 
decomposition of the quark phase within ^ 13 fm/c. The system moves on a path along the 
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phase boundary. In our case we met the situation (S/ A)'^'^^ / {S / A)^'^ > 1 which leads to a 
reheating, i.e. [S/A)'^'-^^ increases with time [p!^, pO| , p^ , to a final value of S/A ^ 140. As a 
consequence the temperature also increases, but only slightly (as can be seen in Fig. |^) . For 
= XQO MeV the system cools, while the specific entropy of the quark drop decreases. 
In this case, the hadronization is incomplete and a strangelet of baryon number A ^ 8 — 9 
is formed. 



Due to the 'strangeness distillery' effect strange and antistrange quarks do not 



hadronize at the same time for a baryon-rich system: Both the hadronic and the quark 
matter phases enter the strange sector, 7^ 0, of the phase diagram immediately. This is 
valid for both scenarios, the higher bag constant (reheating) as well as the lower bag con- 
stant (cooling). The effect can be seen from the high values of the strangeness fraction in 
the quark phase, shown in Fig. ^ for the higher bag constant B^^^ = 235 MeV. At the late 
stage of the evolution, as the strange chemical potential increases, the hadron phase reaches 
positive fs values as well. Fig. ^ shows also the charge to mass ratio of the two phases. 

The rapid expansion of the system leads to changes of the chemical potentials which 
reflect on the differential hadron production rates. Fig. |^ shows the particle rates dNi/dt 
for different hadrons as functions of time for the initial condition A^** = 100, S/A = 45, 
/, = and a bag constant of 5^/^ = 235 MeV (as in Fig. |^). The particle rates decrease in 
general due to the shrinking of the system, since evaporation is proportional to the surface. 
However, the differences in the time dependences of different hadron rates are considerable: 
proton and deuteron rates drop very fast due to the decreasing quark chemical potential, 
while the antiproton rate even increases (for the same reason). K~ and Q production profits 
from the high strange quark chemical potential at the late stage of the evolution. 

Fig. 1^ shows the time evolution of various particle ratios for the same calculation. Most 
of the ratios - the final values of which are given in Figs. ^ - change by factors ^ 2 — 5 in 
the course of the hadronization and evaporation. Therefore, the thermodynamic parameters 
of a certain stage of the evolution, e.g. of the initial stage, cannot be deduced from the final 
integrated particle ratios in the present model. For the bag constant of 5^/^ = 235 MeV the 
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temperature stays rather constant (which is not required by the model). Therefore, ratios 
which do not depend directly on the chemical potentials, like ri/n^ and (p/{p + u), do not 
change significantly. 

Table | shows the absolute abundances of different hadron species for different initial 
conditions of the hadronizing quark-hadron system. We assume an initial baryon number 
of Ab = 208 + 208 (Pb+Pb). However, the final particle ratios do not depend on this, 
because the hadrons yields in this model scale with As- Thus one can easily re- normalize 
the model prediction for any system size, i.e. the number of participant nucleons. Keep 
in mind that isospin symmetric systems are assumed. For B^^^ = 235 MeV the system 
hadronizes completely, while the specific entropy of the quark phase rises from S'/A*"** = 30 
(40, 50) to S/A-^*""' = 109, (130, 153). For B^/^ = 160 MeV a cold strangelet of mass A ^ 35 
(S/A*™* = 40) or A ^ 42 (S/A*™* = 150) emerges. The table shows that antibaryon yields 
are very sensitive to the bag constant, whereas the pions depend strongly on the initial 
specific entropy. 

The K/tt ratio from the dynamic hadronization and the 'favorable' bag constant of 
^1/4 _ MeV turns out to be too high (this is also a problem of the static fit). From 
the preliminary plots in P5| one can read off a K^/h~ ratio of 0.11-0.12 for central Pb+Pb 



collisions {4tt yields), whereas the model yields a value of K^/h^ ~ 0.15 for the three 
choices of initial specific entropy. It is possible to improve this ratio by assuming a lower 
bag constant. For B^^^ = 160 MeV the model renders K^/h~ ratio of ~ 0.1, as can be seen 
in Table |. It can be argued that the numbers of kaons and pions as the most abundant 
particles refiect the actual thermodynamic conditions much more accurately than other 
hadron species. 

Fig. 1^ shows calculated multiplicities of various hypermatter clusters for B^^^ = 235 MeV, 
Ab = 416 and S/A = 40. The penalty factor ~ q^i^-'^)/'^ suppresses the abundances of heavy 
clusters during the hadronization process, which is refiected by the final yields: metastable 
hypermatter can only be produced with a probability p < 10^^ for A> 4 (e.g. a {2S^, 2S°} 
object). 



Can one discriminate the static thermal and the hadronizing source from observations 
of hadron abundances? Obviously, the final particle ratios are nearly identical and therefore 
not suitable for discriminating the two scenarios. The strong change in the time dependent 
particle rates (Fig. H) reflects different average freeze-out times of the particle species. These, 
in turn, correspond to different average freeze-out radii, since the quark drop shrinks — in 
this scenario — during the hadronization process. Both quantities, which characterize the 
size and the lifetime of the particle emitting source, are in principle accessible by means of 
Hanbury-Brown-Twiss-analyses [^. This concept is used extensively for high energy heavy 
ion collisions (for a review see ||2^). 

In Table |I| the average freeze-out times of different hadron species resulting from a 
hadronizing QGP drop (for different initial conditions) are listed. Only the directly produced 
hadrons are taken into account. The true freeze-out times are generally larger due to the 
finite lifetime of the resonances, which contribute to specific hadron yields. However, the two- 
particle correlations within the source — imposed by the quantum statistical momentum 
distributions — are lost through these decays. Thus, the measured correlation function 
does not represent a simple Fourier-like transform of the particle distribution inside the 
emitting source, if the contributions from resonance decays are dominant. As can be seen 
from Table |l| only a fraction of the integrated, final yield of protons, kaons and pions are 
directly emitted from the hadronizing system. A major part stems from the decay of higher 
resonances. This makes the accessibility of the freeze-out radii or times by means of HBT 
analysis questionable. However, the presented results suggest that the observed freeze-out 
time will be smaller for than for K~ , and it will also be smaller for p and A than for 
p and A, respectively. The corresponding radii behave the opposite in this model, because 
the system shrinks. 
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IV. CONCLUSION 



Experimental particle ratios at SPS energies are compatible with the scenario of a static 
thermal source in chemical equilibrium p. A simple dynamic hadronization model, however, 
reproduces the numbers equally well. It yields the following conclusions: a high bag constant 
of B^^^ > 200 is mandatory. These values of B would exclude the existence of stable, 
cold strangelets, as hadronization proceeds without any cooling but with reheating. In 
central collisions of S+Au(W,Pb) a specific entropy per baryon oi S/ A = 35 — 45 is created. 
Due to strangeness distillation the system moves rapidly out of the T,hb plane, into the 
/is-sector. The quark chemical potential drops during the evolution, the strange quark 
chemical potential rises. Final /x^, /z^ values of 1/3 resp. 3 of the values of the static fits are 
reached for t — > t freeze-out- We have presented the model calculations of absolute yields of 
different hadron species for Pb+Pb collisions at SPS energies. The average freeze-out times 
of different hadron species differ strongly, which might be observable via HBT-analyses. 
Strangeness to baryon fractions of ~ 1 — 2 suggest that 'A-droplets' or even 'S^-droplets' 
form the system at the late stage. 
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TABLES 

TABLE I. Hadron abundances for Pb(160GeV/u)Pb (A*^^* = 416, fs = 0) and different initial 
conditions of the hadronizing quark drop. The numbers for the total yields (including feeding from 
higher resonances) are given. 
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TABLE II. Average freeze-out times (in fm/c) of different hadron species for Pb(160GeV/u)Pb 
(^j^intt _ ^-^Q^ pnit _ Q-j fpj^g results for different initial specific entropies of the hadronizing quark 
drop and different bag constants are shown. The fraction of directly produced hadrons divided by 
the total yield (including feeding) is given in brackets. 
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FIG. 1. Particle ratios in the static equilibrium scenario with the parameters (temperature, 
chemical potentials) as indicated. The crosses denote the resulting values if contributions due to 
the decay of higher lying resonances are ignored. The circles include these effects. Data from 
various experiments as compiled in are also shown. 
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FIG. 2. Final particle ratios in the non-equilibrium scenario with initial conditions j4g'* = 100, 
^^^mit _ jmit _ Q g^^^ ^g^g constant B^^^ = 235 MeV. The crosses denote the resulting values 
if contributions due to the decay of higher lying resonances are ignored. The circles include these 
effects. Data from various experiments as compiled in 0] are also shown. 
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FIG. 3. Final particle ratios (including feeding) in the non-equilibrium scenario with initial 
conditions ^g'' = 100, /™* = and bag constant B^^'^ = 160 MeV. The initial specific entropy is 
S/jd™'* = 40 (diamonds) and S'/yl™'* = 150 (squares), respectively. Data from various experiments 
as compiled in |Q are also shown. 
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FIG. 4. Time evolution of the quark chemical potential and the strange quark chemical poten- 
tial for the parameter set of Fig. |2| (here and in the following pictures we start the hadronization 
at an initial time of to = 5 fm/c). The diamond denotes the chemical potentials extracted with 
the static fit (Fig. |). 
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FIG. 5. Hadronization path of the system (initial conditions as in Fig. in the projected T—fiq 
plane. The numbers denote the time in fm/c. Shown is also the phase boundary of a QGP and the 
hadron gas for i?^/^ = 235 MeV and the path of constant specific entropy in a pure quark phase 
with S/A = 45. The full circle indicates the Hq and T values of the static fit. 
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FIG. 7. Particle rates as functions of time (initial conditions as in Fig. §). 
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FIG. 8. Time evolution of various particle ratios (including feeding) for the parameter set of 
Fig. |2|. Here the ratios of the particle rates at a certain time t are depicted. 
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FIG. 9. Calculated multiplicities of hypermatter clusters from a hadronizing QGP with 
j^nit ^ 415^ = 40, /i™* = and a bag constant of B^/^ = 235 MeV: 

a) (m = 2020 MeV), b) {E-,E^}, c) ^He, d) {4A}, e) {25", 25°}, f) iHe, g) j^lHe, 
h) {2n, 2A, 2S-}, i) {2A, 22°, 22"}, j) ^0^^^, k) A = 8, S = 0, I) A = 8, S = -4, 
m) A = 8, S = -8,n) A = 8, S = -12, o) ^ = 8, 5 = -16 
{—S gives the number of strange quarks). 



24 



